Perspective on rainbow-ladder truncation 
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Prima facie the systematic implementation of corrections to the rainbow-ladder truncation 
of QCD's Dyson-Schwinger equations will uniformly reduce in magnitude those calculated mass- 
dimensioned results for pseudoscalar and vector meson properties that are not tightly constrained 
by symmetries. The aim and interpretation of studies employing rainbow-ladder truncation are 
reconsidered in this light. 
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That the lightest bound states supported by QCD are 
members of a J PC = h isotriplet owes to the near 
chiral symmetry of the light-quark Lagrangian and the 
pattern by which it is dynamically broken. The particu- 
lar behaviour of the pion mass; viz., m 2 oc rh, where rh is 
the renormalisation-group-invariant light-quark current- 
mass, and numerous other fascinating properties can be 
derived in QCD through a careful consideration of the 
axial-vector Ward-Takahashi identity and the gap equa- 
tion; e.g., Ref. [3, H, H, 0> Hi- In relation to the lat- 
ter, Dyson-Schwinger equation studies [1,0] and simula- 
tions of the lattice-regularised theory Q have established 
conclusively that dynamical chiral symmetry breaking 
(DCSB) is a reality in QCD. These facts have been 
used to provide a satisfactory understanding of the pion's 
mass. 

Such is not the case, e.g., for the pion's radius nor its 
weak decay constant. They are both well known experi- 
mentally [H: TV = 0.672±0.008fm; f w = 130.7±0.4 GeV, 
but the problem of understanding them within QCD re- 
mains. Their reliable calculation requires more detailed 
knowledge of strong interaction dynamics than does m v . 
Furthermore, one must employ a formulation of the quan- 
tum field theory bound state problem that veraciously in- 
corporates DCSB, and its corollaries, and simultaneously 
preserves the vector Ward-Takahashi identity. 

This can be achieved in a Poincare covariant and sym- 
metry preserving treatment of quark-antiquark bound 
states through the Bethe-Salpeter equation (BSE) 



T tu (k-P)= f 

Jq 



(1) 



where: k is the relative and P the total momentum 
of the constituents; r,. . . ,u represent colour, Dirac and 
flavour indices; x(<l>P) '■— S(q + )T(q; P)S(q ), with 
T(q; P) a given meson's Bethe-Salpeter amplitude and 

q± = q±P/2; and J represents a translationally invari- 
ant regularisation of the integral, with A the regularisa- 
tion mass-scale Q, E3- In Eq. ©, S is the renormaliscd 
dressed-quark propagator and K is the fully amputated 
dressed-quark-antiquark scattering kernel. The product 
SSK is a renormalisation point invariant. 

The dressed-quark propagator appearing in the BSE's 



kernel is determined by the renormalised gap equation 1 



Sip)- 1 = Z 2 (z 7 .p + m bm ) + £(p), 



(2) 



S(p) = Z 1 J g 2 D^(p-q)— 7fl S(q)T a u (q,p),(3) 

where D M „(fc) is the dressed-gluon propagator, T u (q,p) 
is the dressed-quark-gluon vertex, and m bm is the A- 
dependent current-quark bare mass. The quark-gluon- 
vertex and quark wave function renormalisation con- 
stants, Zi i{C , -A 2 ), depend on the renormalisation point, 
(, the regularisation mass-scale and the gauge parameter. 
The gap equation's solution has the form 



s(pY 



l7 -M( P 2 ,c 2 ) + iV,C 2 ). 



(4) 



and the mass function M{p 2 ) = B(p 2 X 2 )/A(p 2 ,( 2 ) is 
renormalisation point independent. The propagator is 
obtained from Eq. ^ augmented by the condition 



s( P y 



<2 



i7 • p + m(C) : 



(5) 



where m(Q is the renormalised mass: Za(Ci A ) m (A) 
= Zi(C, 2 , A 2 ) m(£), with Z± the Lagrangian mass renor- 
malisation constant. The chiral limit means rh = 0. 

The Bethe-Salpeter and gap equations are Dyson- 
Schwinger equations (DSEs), for which a nonperturba- 
tive Poincare covariant and symmetry preserving trunca- 
tion scheme is available [ill, El • The leading-order term 
is the so-called rainbow-ladder truncation. We assume 
it is an accurate tool for calculating meson properties 
in the heavy-quark limit because for relevant momenta 
the contribution from additional diagrams diminishes as 
the curre nt-q uark mass of a bound state's constituents 
increases [13|. Furthermore, in equal-mass-constituent 
pseudoscalar and vector meson channels it can be seen al- 
gebraically that contributions from subleading diagrams 
interfere destructively amongst themselves owing to the 
axial- vector Ward-Takahashi identity [H, EE 13 • 



1 In our Euclidean metric: {'y^iOv} — 25^^: 
74717273; a ■ b = J2i=i a i b i' and p n timelike 



=>■ P 2 < 0. 
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These facts have successfully been exploited through 
the application of a renormalisation-group-improved 
rainbow-ladder DSE model [l(| EH] to a wide variety of 
phenomena involving it- and p-mesons (l6j . The heart of 
the model is an Ansatz for the Bethe-Salpeter kernel: 



K t r u s (q,k;P) = -G((k-q) 2 ) 
X D^(k q) 
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(6) 



wherein D^^ e (£) is the free gauge boson propagator and 

nana El 
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Z 2 s 



■■ C{oj, rh) 
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s/A 2 



QCD 



e -s/k 2 A?] 

^( S ) , (7) 



with -F(s) 



[1 - cxp(- s /A 2 )]/ s , A t 



A. 



(4) 



= 1.0 GeV, r 
0.234 GeV. 



e' - l, 7m = 12/25 and Aqcd -- Us 

This form expresses the interaction as a sum of two 
terms. The second ensures that perturbative behaviour 
is correctly realised at short range; namely, as written, 
for (k — q) 2 ~ k 2 ~ q 2 > 1 — 2 GeV 2 , if is precisely as 
prescribed by QCD. On the other hand, the first term 
in G(k 2 ) is a model for the long-range behaviour of the 
interaction. It is a finite width representation of the form 
discussed in Refs. [201 ]. 

Given a truncation of the BSE's kernel there is a unique 
gap equation that ensures the Ward-Takahashi identities 
are automatically satisfied [HI E2, 0) B3 • In connection 
with Eq. (|BJ) that is a rainbow gap equation; viz., Eq. ((2J) 
with 



g((p-q) 2 )D^(p- 



<l)^l^S{q) — -i u . (8) 



Having selected a DSE truncation and specified the 
behaviour of the interaction, one can straightforwardly 
proceed to calculate hadron properties. There are ap- 
parently two active parametric quantities in Eq. ((7|); viz., 
C and u), which together determine the integrated in- 
frared strength of the rainbow-ladder kernel. They were 
fitted to a selection of physical observables in Ref. [15| ]. 
However, in reconsidering that fit it was noted [18j that 
the predictions are approximately independent of ui on 
the domain uj e [0.3,0.5]. Hence, there is truly only 
one parameter and experimental values are satisfactorily 
reproduced with C{u>,m) — const. = 0.37. This obser- 
vation has since been much exploited. Nevertheless, we 
choose to reconsider this approach and explore whether 
there is a more appropriate way to determine C{u>,rh). 

To develop an answer we remark that in pseudoscalar 
and vector channels the net effect of corrections to the 
rainbow-ladder truncation is attraction. Hence at the 
physical light-quark current-mass the truncation should 
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FIG. 1: Solid line: expected evolution in rainbow-ladder trun- 
cation of the p-meson mass with pion mass-squared, Eq. © . 
Results from simulations of lattice-regularised QCD are also 
depicted: ETMC [H; RBC/UKQCD and MILC [Uj], 

plus the results in Ref. [26[ with an analysis and chiral extrap- 
olation [27;]. In studying this figure it is important to bear in 
mind that Eq. is a calculated overestimate of m p (in^). 



yield a value for m p that is larger than experiment. Ex- 
tant studies indicate that a reasonable value is mJJ := 
m^ L (m = 0) = 0.99 GeV because the resummation of 
nonrcsonant diagrams 0, 0, H3 produces ~ 100 MeV 
of attraction and pseudoscalar meson loops provide a 
further ~ 120 MeV, in addition to generating a width 
[22] ]. Corrections to rainbow-ladder vanish with increas- 
ing current-quark mass. For example, at rh such that the 
mass of a pion-like pseudoscalar is 0.63 GeV, corrections 
to rrip from meson loops are already negligible. Further- 
more, nonresonant contributions are 70% of their chiral- 
limit magnitude; at m, = 0.9 GeV that drops to 50%; 
and they vanish completely in the heavy-quark limit. We 
express this behaviour through (see Fig.[T]): 



m 2 (m 2 ) 



= 1 



(rn^/m ) 4 



0.6 + {rriTr/mPp) 



(9) 



Since vector meson masses are observable, well known 
and not constrained by symmetries, we determine 
C(w,7n) by requiring that, with the interaction model 
defined in Eq. (|7|) , rainbow-ladder truncation reproduce 
Eq. ([9]). The procedure is straightforward. For a given 
current-quark mass, and w in Eq. ([7]), the rainbow-ladder 
truncation produces values for m p and via Eqs. |T]), 
@, ^ and ||5J). C(u),m) is chosen such that these val- 
ues are related through Eq. ([9|) . The result is depicted 
in Fig.Elwhich also illustrates clearly the observation 
of Ref. 18]; viz., for each rh there is a domain of uj on 
which Eq. ([9]) is preserved with C(u), rh) w constant. The 
shaded band in Fig. [2] highlights that domain. We will 
subsequently take this band to define a domain of in- 
ternally consistent interaction tensions [28| and report 
physical quantities for all values of ui within this domain. 
The variation in results should be taken as an indication 
of the model-dependent uncertainty in our predictions. 
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FIG. 2: C(w,m) in Eq. © determined by Eq. Q. Five 
current-quark masses are shown explicitly. The shaded re- 
gion demarcates a domain of minimal sensitivity to this 
model parameter. The dashed almost-vertical curve indicates 
a(m) = 0.38 + 0.17/[1 + tfi/jfto], rh = 0.12 GeV. 



An accurate interpolation within this band is 



Ci(cj - u(m)) „ 
C{uj,m) = C Q + - ; - - ; ^ 2 „ 2 , a; = m/m , 



1 + C 2 x + Cf x 
Ci(£) = 0.86(1 - 0.15£ + (1.50f) 2 



(10) 



(2.95f) 3 ). 



with L>(m) the m-dependent midpoint of the shaded do- 
main in Fig.[2]and C = 0.11, C 2 = 0.885, C 3 = 0.474. 

Figure [3] illustrates the impact of this analysis on the 
behaviour of the functions in Eq. ^ , which characterise 
the dressed-quark propagator. For relevant momenta, 
corrections to the rainbow truncated gap equation are 
suppressed by increasing to. Furthermore, while they are 
always negligible for p 2 > 3 GeV 2 , corrections suppress 
M(p 2 ) and Z(p 2 ) for p 2 < lGeV 2 , while maintaining 
their support at intermediate momenta [HI]. 

We have now determined the dressed-quark propagator 
and meson Bethe-Salpeter amplitudes, and their current- 
quark mass dependence. Physical quantities can there- 
fore be calculated; e.g., the vacuum condensate [l|, 



(qq)° c = hm Z 4 (( 2 , A 2 ) N c tr D 



A^oo 



S°(q,C), (11) 



where the superscript "0" denotes a quantity calculated 
in the chiral limit, and the it and p decay constants [H.[29|: 




FIG. 3: Top panel - mass function and bottom panel - wave 
function renormalisation, Z(jp,Q = l/A(p,Q, calculated with 
C = 2.9 GeV and rh such that m n = 0.63(1) GeV. For M(p 2 ), 
the lower (upper) edge of the bands corresponds to the largest 
(smallest) u) value in the domain demarcated in Fig. [2]- this 
pattern is reversed for Z(p, £). Circles in each panel are results 
from simulations of lattice-regularised QCD 



TABLE I: Row 1: Results calculated with rh = 6.1 MeV, for 
which — 0.138 GeV. Error describes sensitivity to varia- 
tion in lu around Q in Fig. [5] Row 2: Experimental [i| or phe- 
nomenological values. All quantities in GeV, except 7V in fm. 
The condensate is quoted at 1 GeV after one-loop evolution 
from the renormalisation scale C," = 19 GeV, as described in 
Ref. [l(j. The axial- vector Ward-Takahashi identity ensures 
that subleading corrections to the pion's mass are < 5% [l3| ]. 
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0.319(2) 
0.236 


0.178 
0.131 


0.283(9) 
0.216 


0.99 
0.77 


0.485(3) 
0.672 



/ W+ P„ = Z 2 tr CD / 757^ (q;P), (12) 



f P +m p = iz 2 tr CD / ^ (q; P) ■ (13) 



In Table U we present results obtained in rainbow- 
ladder truncation using Eqs. ([7]) and (|10[) at a current- 
quark mass for which to^ = 0.138 GeV. Expressed in 



units of mass, the tabulated quantities overestimate ex- 
periment by 34 ± 4 %. This uniform response to a delib- 
erate increase in the interaction tension is internally con- 
sistent; viz., a single cause produces a uniform effect and 
that permits a single remedy. NB. Attempts to estimate 
the contribution from pseudoscalar meson loops alone 
produced uncorrected values of [— (qq)®]^ ~ 0.28 GeV 
and ~ 0.17 GeV [30l [3l|. comparable to Row 1. 
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In Fig. [4] we depict the calculated current-quark mass 
dependence of the pion's leptonic decay constant. A ma- 
terial similarity with Fig.Q] is apparent. NB. The dis- 
cussion associated with Eq. © should be borne in mind 
when contemplating this figure. 

In rainbow-ladder truncation the electromagnetic pion 
form factor is determined from 13211 



(Pi + P2)^(Q 2 ) = tr CD J f „+(k; -pa) S(k+) 

x r 7I+ (k+; Pl )S(kl)tT^k+,kZ)S(kZ), (14) 

where fcf = k + ap 1 /2 + (3Q/2 and p 2 = Pi + Q. The 
only quantity in Eq. (|14[) that we have not yet encoun- 
tered is the dressed-quark photon vertex, r p (£i,£ 2 )- It 
is obtained from an inhomogeneous Bethe-Salpeter equa- 
tion and, in a symmetry preserving DSE truncation, the 
solution satisfies the vector Ward-Takahashi identity. 

The definitive rainbow-ladder calculation of Ft T {Q 2 ) 
solved the vertex equation and employed the result in 
Eq. (|14|) [331 ] . It can be argued from this study that a re- 
liable approximation to the result for the pion's charge ra- 
dius in rainbow-ladder truncation is provided by [33, H3] 



' ttBC 1 2 



6 g P w(fi) 



9p 



(15) 



The first term in Eq. (|15|) . r^sc, is the nonresonant 
contribution to the radius. It is obtained from Eq. (TT4"]) 
with the Ball-Chiu Ansatz for 1^(4,4) HH- The 
second term represents the p-meson contribution to the 
radius. That is explicitly excluded from Eq. (| 14[) when 
one employs the Ball-Chiu vertex calculated from the 
solution of the rainbow gap equation. The quantity 
g p = \f2m p / f p is the p° — > e + e _ coupling. The function 
gpmriQ 2 ) describes the pirn vertex, which in rainbow- 
ladder truncation is given by a straightforward analogue 
of Eq. (Ti"4"|) . Experimentally, g p7r7r (Q 2 = — mi) = 6.14 
but Eq. (fTS]) requires an off-shell-/? coupling, g P T T - K {Q 2 — 
0). We calculate that with the on-shell p- meson Bethe- 
Salpeter amplitude, modified such that each sublead- 
ing Dirac structure acquires a suppression factor: [2 + 
Q 2 /m 2 }^ 1 / 2 . The factor e _e in Eq. (fTS")) accounts for 
any additional effect owing to an off-shell p-meson. 2 Its 
value is determined by requiring that Eq. (|15p repro- 
duce the DSE curve in Fig. 6 of Ref. [36[ , with the re- 
sult q — 0.2m T . This analysis gives g P 7r7r(0)e~ e w 
constant = 2.75(40) on the relevant (w,m,) domain: 
larger u>, smaller constant value, and vice-versa. 

Thus specified, Eq. (TTSJ) yields the charge radius in Ta- 
ble[T] and the evolution with m T depicted in Fig. [5] On 
the material (w,m ff ) domain the nonresonant contribu- 
tion provides approximately 50% of the reported value 
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FIG. 4: Shaded band: calculated mass dependence of the 
pion's leptonic decay constant. The band delimits the range 
of results possible as uj varies over the domain in Fig. [2] Re- 
sults from simulations of lattice-regularised QCD are also 
presented: ETMC HI; RBC/UKQCD 0]; HPQCD H; 
NPLQCD [H; MILCgcj; QCDSF 0; and CERN [HEl, 
a compilation drawn from Ref. 




m n [GeV] 

FIG. 5: Shaded band: calculated pion-mass dependence 
of the pion's mean-square charge radius. The ^-variation 
band is narrow owing to cancellations between effects on the 
two terms in Eq. (|15[1 . Results from simulations of lattice- 
regularised QCD are also presented: quenched, triangles and 
squares 45, 4(|; and unquenched, circles [47| . 



of r 2 . Furthermore, on this domain we calculate r^f^ = 
0.44(1); viz., a constant independent of the current-quark 
mass. The experimental value is 0.445 ± 0.007. As could 
be anticipated from Ref. [43], with rsc alone the calcu- 
lated value is only ~ 0.32. It is noteworthy that an es- 
timate of the contribution to the charge radius owing to 



pseudoscalar meson loops gives [3j 



7T loop 



0.07(2) fm 2 



2 A single-pole vector meson dominance Ansatz assumes 
9pmr{Q 2 = 0) = g p , Q = 0. The latter is also assumed in Ref. [33j. 



at to w = 0.14 GeV. This must simply be added to our re- 
sult because the p-meson term in Eq. (|15|) has zero width. 
The result is r„ = 0.55(2). Naturally, as m n increases the 
7r-loop contribution diminishes rapidly in importance. 
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We have explained that in connection with light-quark 
systems, and those of the physical qualities of the pseu- 
doscalar and vector meson bound states they constitute 
which arc not tightly constrained by symmetries, the 
rainbow-ladder truncation of QCD's DSEs should pro- 
duce results that, when measured in units of mass, are 
uniformly « 35% too large. The systematic implemen- 
tation of corrections will then shift calculated results so 
that reliable predictions and agreement with experiment 
can subsequently be expected. In this way one can arrive 
at a veracious understanding of light-quark observables. 

Our study also explains why, when employed as a 
means of modelling QCD, it is possible to tune the para- 
metric elements in a rainbow-ladder truncation such that 
a wide range of light-quark pseudoscalar- and vector- 
meson observables can successfully be correlated 



That rainbow-ladder is the first term in a systematic 
truncation guarantees a uniform response to corrections 
and hence allows for their effects in pseudoscalar and 
vector channels to be expressed in large part through 
simple parameter modifications. Naturally, however, fea- 
tures such as decay widths and associated nonanalytic- 
ities, which are essentially tied to hadron loops, cannot 
be realised in this way. 
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